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C 1 Abstract 



Chaotic regimes of the microwave energy absorption are experimentally 
. observed and analyzed for two-dimensional metallorganic antiferromagnet 

£^ . (NH3)2(CH2)4MnCi4 at low temperatures under the conditions of nonlinear 

antiferromagnetic resonance. Relaxation oscillations of energy absorption are 
investigated in detail. Their frequency spectra, frequency-amplitude char- 
acteristics, and dependences of absorbed power on driving power and static 
, magnetic field are studied. It is shown that the dynamics of relaxation oscil- 



X 



lations undergoes a transition to chaos by "irregular periods". Peculiarities 
of the transition are described consistently. Among other things, the condi- 
tions for the emergence of energy absorption regimes with a spike-like and a 
saw-tooth signal structure are determined, and the characteristics of chaotic 
oscillations such as the dimensions of strange attractors are calculated. The 
chaotic dynamics is found to be high-dimensional with a large contribution 
from noise which is of deterministic origin in the antiferromagnet under in- 
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vestigation. 
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INTRODUCTION 



Chaotic resonant phenomena in magnets have become an object of intense experimental 
studies in the last decade. These investigations were stimulated by the progress made in 
the mathematical theory of chaos predicting the universal character of chaotic phenomena 
irrespective of the nature of the physical object being studied and demonstrating a nonlinear 
behavior .i^il 

Magnetic compounds possessing the properties required for the emergence of nonlinear 
oscillations include first of all the crystals exhibiting an extremely weak relaxation of spin 
excitations. Yttrium-iron garnet (YIG) with a low threshold for a parametric excitation of 
spin waves even at room temperature has been studied most thoroug hiy.!H!£K! Since YIG 
behaves as a ferromagnet in the magnetic respect, nonlinear chaotic effects were studied, as 
a rule, under the conditions of ferromagnetic resonance (FMR) in a transverse as well as 
longitudinal driving fields. 

The total number of investigated nonlinear magnet is not large, and some of them exhibit 
nonlinear properties only at low temperatures of the order of a few kelvins, at which phonons 
are frozen out, and their interactions with magnons becomes very weak. 

The effective dimensionality of crystals plays an important role for the suppression of 
relaxation processes. Stepanov et alJlUll investigated the class of metallorganic compounds 
that are quasi- two-dimensional ferro- and aniferromagnets in the magnetic respect. It was 
found that low- dimensional magnets at low temperatures go over to a state with an anoma- 
lously low spin-lattice relaxation virtually with a threshold, which makes it possible to excite 
spin waves parametrically at microwave pumping power of the order of a few milliwatts!! 
Among other things, it was found that in addition to YIG, nonlinear ferromagnetic crystals 
include metallorganic compounds with a structure similar to the (CHaNHa^CuC^ crystal.ll! 
On the other hand, it was established that chaotic oscillations are generated in the crys- 
tals CuCl 2 2H 2 0,! CsMnF 3 ,l§ (CHaNHs^CuClJ under the conditions of antiferromagnetic 
resonance (AFMR). 
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The range of nonlinear effects that have been discovered and thoroughly investigated in 
ferro- and antiferromagnets is quite large. These include spin-wave instabilities (Suhl insta- 
bilities of the first- and second order) ,0 auto-oscillations of absorbed microwave power ,@H 
and the observation of three known scenarios of a transition to chaos: by period doubling 
(Feigenbaum scenario), quasiperiodicity, and intermittencyB^ii 

Apart from the interpretation of these nonlinear effects and the determination of the 
conditions for their observation, it was found that real magnetic crystals can demonstrate 
a more complex pattern of transition to chaotic regimes in resonance experiments. None 
of the known scenarios is realized in pure form in such cases,@ and we must consider new 
mechanisms of chaotizationBiil 

Hartwik et al.il were the first to discover long ago the so-called relaxation chaotic oscil- 
lations of microwave power in YIG with which a new scenario of a transition to chaos by 
irregular periods has been associated in last decade .i'iffi This effect lies in the emergence, 
instead of purely periodic auto-oscillations, of irregular chaotic bursts of absorbed power in 
the form of spike-like peaks or pulses with a steep leading front and relaxing rear front under 
certain conditions of magnetic resonance upon an increase in the pumping power. Theo- 
retical approaches to the description of such oscillations and mechanisms of their formation 
were made in Refs. |9 p9[ , [30i but a systematic analysis of temporal series of experimental 
signals as well as of the results of numerical simulation of relaxation oscillations was carried 
out only recently.0H 

In this work, we study experimentally the regimes of chaotic behavior of the microwave 
power absorbed in a two-dimensional easy-axis antiferromagnet (NHs^CH^^MnCU under 
the AFMR conditions. This compound is a typical representative of the family of layered 
Heisenberg antiferromagnets [NH3-(CH 2 ) m -NH 3 ]MnCl4 (2CmMn) studied by Stepanov et 
The structure of these metallorganic crystals is formed by almost quadratic layers of 
magnetic ions in the octahedral environment of chlorine ions between which long chains of 
alkylene-ammonia molecules are located. The small value of interlayer exchange associated 
with a large separation between the spins of adjacent layers leads to a quasi- two-dimensional 



behavior of these systems. At temperature T/v = 42.6 K,E3 the compound 2C4Mn is trans- 
formed into the antiferromagnetic states with the easy magnetization axis directed at right 
angles to the planes of the layers. A detailed analysis of linear antiferromagnetic resonance in 
2C4Mn revealedS that this compound has a four-sublattice noncollinear antiferromagnetic 
structure with a weak ferromagnetic moment. The antiferromagnetism vector of each layer 
is deflected successively from the normal to the plane of a layer through an angle of ±16° so 
that the total vector is perpendicular to the layer, and the weak antiferromagnetism vector 
lies in the layer. According to estimates, the strength of interaction between the layers 
(26 Oe) is extremely small as compared to both the intralayer exchange (2H e ~ 1360 kOe) 
and the intralayer uniaxial anisotropy (H a ~ 0.8 kOe) so that 2C4Mn can be regarded as an 
almost two-dimensional antiferromagnet. 

At low temperatures (of the order of a few kelvins), the related compound 
(NH3C2H5)2MnCi4 (lC2Mn) revealed a number of interesting nonlinear phenomena in the 
behavior of the absorbed microwave power, e.g., the emergence of auto-oscillations and 
chaos.0 Here we carry out a systematic analysis of chaotic regimes of antiferromagnetic 
resonance in a 2C4Mn crystal. For a driving power below 5 mW at a temperature below 
2.18 K, we observed a nonlinear absorption of the microwave field and the emergence of 
relaxation oscillations with typical (extremely low) average frequencies of the order of a few 
hertz. These oscillations were recorded and analyzed as temporal series of data with the 
help of an analog- digital device and computer programs, which made it possible to describe 
in detail description of the scenario of a transition to chaos by irregular periods. As a re- 
sult, we have analyzed qualitative changes in the behavior of temporal series of absorbed 
power as a function of variation of the parameters of static and varying magnetic fields and 
carried out the Fourier analysis, obtained the spectra of oscillations, studied the structure 
of strange attractors of chaotic regimes, and calculated the quantities characterizing chaotic 
dynamics, i.e., dependences of frequencies of auto-oscillations oscillations on the driving 
power (in particular, we determined their period doubling threshold) and the dimension of 
chaotic attractors, and discussed the origin and role of noise in relaxation oscillations as 
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well as a possible theoretical model for describing relaxation oscillations in two-dimensional 
ant if erromagnet s . 



EXPERIMENTAL TECHNIQUE 



Single crystals of (NH 3 ) 2 (CH 2 )4MnCl4 were grown at room temperature from a saturated 
aqueous solution in the form of thin rectangular plates with clearly manifested lateral faces 
and with a typical size 5 x 5 x 0.3 mm. The compound (NH 3 ) 2 (CH 2 )4MnCl4 possesses a 
monoclinic symmetry of crystal lattice with the space group P2i/b^ Organic chains of 
NH3(CH 2 )4NH3 separate two-dimensional almost square layers of octahedra Mn-Cl6- The 
unit cell parameters are a = 10.77 A, b — 7.177 A, c = 7.307 A. Experiments were carried 
out on a reflection spectrometer with a pumping frequency 70.39 GHz at a temperature 
below 2.18 K. We used a cylindrical resonator with the Q- factor ~ 1000. The sample was 
placed in the resonator region with predominant parallel polarization of external static and 
rf fields H \\ h. 

In resonance experiments, the field is usually applies along the easy axis of the crystal. 
With such an orientation, the splitting of AFMR branches follows the law 



For the frequency mentioned above, the resonance conditions are satisfied for the lower 
frequency branch tu_, which was observed in our experiments. The maximum power of the 
source was 5 mW. The applied magnetic field was scanned along the contour of the AFMR 
line, and the driving power was varied from to —20 dB. The magnetic field orientation 
relative to the anisotropy axis and equilibrium directions of antiferromagnetism vectors in 
adjacent planes also varied. It was found that the most intense absorption corresponds to 
the symmetric orientation of the field along the crystallographic axis b. 

In all experiments, low-frequency modulation of electromagnetic field of frequency 50 Hz 
was observed. This frequency had to play the role of the reference frequency in our experi- 
ments. It was found later that these oscillations participated in all nonlinear processes, and 




(1) 
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the emergence of their higher harmonics was regarded as a natural criterion of the emergence 
of nonlinearity in resonance effects. The reflected signal after detection in an analog- digital 
device PC ADDA- 14 with a 14- bite resolution was transformed into a computer data file. 
These temporal series were subsequently analyzed by using the standard and original packets 
of programs created for a quantitative analysis of chaotic phenomena. 

DISCUSSION OF MAIN RESULTS 

A typical form of resonant curves for the antiferromagnet 2C4Mn at T = 1.8 K are 
shown in Fig. |l[ For a low (less than — 15 dB) microwave field power, a typical pattern 
from a linear AFMR is observed, i.e., two lines from two centers (neighboring planes). The 
separation between the peaks on the resonant curve can vary depending on the orientation 
of the static magnetic field, and the lines can coincide when the field is directed along 
the crystallographic axis b. For a power exceeding — 15 dB, free relaxation oscillations are 
generated in the range of external magnetic fields near the peak of the high-field line above 
as well as below the resonant field H = 9.34 kOe (this region is shown by the bold line on the 
upper curve in Fig. 1). The amplitude of these oscillations increases with power until they 
become chaotic. We shall consider this regime in detail later, and now we pay attention to 
another effect associated with instability of resonance at high pumping levels. As the driving 
power increase above — 5dB (see two lower curves in Fig. 1), a jump and a discontinuity of 
both resonant lines are observed with a considerable hysteresis in the static magnetic field 
as we move towards higher and lower fields respectively. 

This phenomenon is well known in the theory of nonlinear resonance and is associated 
with the dependence of the frequency of nonlinear oscillations on their amplitude. In mag- 
nets, this effect is manifested in that the resonant curve must become asymmetric and 
multiple- valued for a pumping field h exceeding the critical value, i.e., the peak must be in- 
clined towards lower or higher fields depending on the type of interaction between magnons. 
In actual experiments, instability is observed upon a change in the static magnetic field, and 
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the resonant curve experiences a discontinuity or a jump. This effect was observed for the 
first time in disks of yttrium-iron garnet single crystals by Weiss.H 

It follows from Fig. 1 that the jump is observed in strong fields in an increasing field, 
while discontinuity takes place in weaker fields in a decreasing field. With increasing power, 
the hysteresis loop increases, and the steepness of lines decreases (the scales of conditional 
units for absorbed power in Fig. 1 are different for the three resonant curves: it decreases 
with increasing amplitude of pumping). 

It was noted above that oscillations of observed power on the segment of the resonant 
curve near 9.34 kOe appear even at very low driving powers of the order of — 15 dB. The 
criterion of a transition to the nonlinear regime is the emergence of the second harmonic 
peak in reference oscillations with frequency 100 Hz. At the point of maximum on the 
resonant curve, this peak exceeds the background noise for a power P > — 15 dB, and first 
spike-like peaks of absorbed power appear at the same instant. 

In the case of a resonant curve with spaced peaks, an increase in the driving power induces 
relaxation oscillations in the vicinity of the second peak also. As the peaks converge, the 
mutual effect of the centers increases, which is noticeably reflected in the form of oscillations 
of absorbed power. In the cases of closely spaced peaks for the value of the field H m = 
8.37 kOe corresponding to a local minimum at the center of the resonant curve, relaxation 
oscillations become irregular even for a low driving power. As the power increases, the 
oscillations become more and more chaotic. The time dependence of the signal typical of 
the entire series of these measurements and its spectrum for the maximum value of power 
are shown in Fig. 0a. 

In order to find out whether such a dependence is a consequence of additive or dynamic 
noise, stochastic process, or is due to a determinate chaos, we varied in the experiments the 
orientation of magnetic field and its magnitude. The variations affected strongly the type 
of oscillations. 

It was found that oscillations become less chaotic for the minimum deviation of the 
magnitude of magnetic field from the extremal value H m . By way of an example, Fig. |2|& 
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shows the time behavior and spectrum of oscillations for H = 8.4 kOe and the maximum 
driving power. 

It also turned out that the degree of chaotization of a signal decreases considerably, and 
its shape changes qualitatively when the magnetic field is directed along the crystallographic 
axis b, when the resonant lines from two centers coincide. In this case, the characteristic 
pattern of the emergence and transformation of relaxation oscillations upon a change in the 
driving power is of the form shown in Fig. [| for H = 8.4 kOe (small deviation of the field 
from the resonant value) and the power P varying from —10 to dB. First spikes of absorbed 
power appear against the background of almost linear oscillations of frequency 50 Hz. For 
small pumping amplitudes, the frequency corresponding to the emergence of spike-like peaks 
is low, and the intervals between them are quite large and vary with an obvious periodicity. 
For a driving power of the order of —6 dB, the signal has the form of a periodic structure of 
spike-like closely spaced peaks. As the pumping amplitude increases further to power values 
of the order of —3 dB, the frequency corresponding to the emergence of spike-like peaks 
changes insignificantly, and subsequently decreases rapidly and becomes virtually equal to 
half the previous value. In this region, the shape of the signal changes qualitatively from 
the spike- like to the saw-tooth, i.e., the change in the regime of chaotic oscillations takes 
place. Figure |] shows the dependence of the fundamental frequency of these oscillations on 
the driving power in the range from -6 to dB (solid circles). The doubling of the period 
of relaxation oscillations can be seen clearly in the figure. In order to plot this dependence, 
we analyzed the spectra of oscillations for fixed pumping levels. It should be noted that 
doubling of this period does not indicate the emergence of subharmonics of the fundamental 
frequency as is usually the case in the Feigenbaum scenario, and corresponds to the change 
from one oscillatory mode to another mode, their fundamental frequencies differing by a 
factor of two. It was proved that this effect is preserved for other values of magnetic field 
which naturally affects the values of frequencies themselves. It should also be noted that 
apart from the main peaks and multiple harmonics, all spectra contain a large contribution 
from noise responsible for the "grass-like" continuous spectrum. In the subsequent analysis, 
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we shall analyze in detail the dynamic and spectral structure of these oscillations and the 
origin of their stochastic form. 

When the power changes in the opposite direction, i.e., the amplitude of pumping de- 
creases (open circles in Fig. the frequency-amplitude dependence exhibits a hysteresis 
with a displacement of the region of period doubling towards lower powers (saw-tooth pulses 
exist up to —4.5 dB). The existence of essentially chaotic modes near a certain fixed values of 
power, in particular upon an increase in the driving power for P = — 1, —2.25, and —2.75 dB 
is an interesting feature of the observed transient process. 

For this reason, it was natural to analyze oscillations for these selected pumping levels, 
but in a wide range of applied magnetic field near the resonance point. We chose the 
pumping level of —1 dB and studied the variation of the shape of the absorbed power signal 
and its spectrum upon a change in the static magnetic field within a few ten oersted near 
the resonant value H r = 8.37 kOe. The direction of the field was maintained along the 
crystallographic axis b. 

The corresponding results are presented in Fig. [5]. It should be noted that relaxation 
oscillations occur against the background of a considerable average absorbed power that 
must make a contribution to the frequency spectrum in the form of a large central peak 
at zero frequency. In all calculations of the spectra analyzed here, this average value was 
subtracted, and hence the given huge contribution to the central peak is absent, which allows 
us to see the detailed structure of relaxation oscillations proper. It should also be noted that 
frequency spectra are given in the form of frequency dependences of the amplitude of the 
Fourier transform of the signal and not as logarithmic spectra of power in order to improve 
detailization. 

Another feature in common of all the spectra considered below is the presence of oscil- 
lations of frequency uq = 50 Hz. These low-frequency oscillations were present as a source 
of reference frequency, but they became involved in free oscillations in view of the nonlin- 
earity of the medium. This follows from the presence of second harmonic with frequency 
v = 100 Hz and the peaks that are algebraic sums of frequencies of fundamental harmonics 
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of relaxation oscillations and the reference frequency. 

Far away from the resonant field, the absorbed power is virtually constant if we disregard 
extremely low background noise in which, however, oscillations with frequency i>q = 50 Hz 
were always manifested (in the frequency spectrum) in our measurements. As the field 
approaches the resonant level, these small-amplitude oscillations become weakly nonlinear 
(second harmonic appears in the spectrum), and nearly periodic spike-like peaks of ab- 
sorbed power corresponding to peaks of the order of a few hertz in the frequency spectrum 
and clearly distinguishable against the "grass-like" background noise appear almost simul- 
taneously. 

A typical example of such a behavior of absorbed power is shown in Fig. |5|a for the field 
value H = 8.51 kOe. It can be seen that periodic relaxation oscillations with a spike-like 
structure have been formed completely. Small anharmonic modulation of peak amplitudes 
is manifested in the frequency spectrum in the form of higher harmonics of the fundamental 
frequency. All the remaining peaks can be identified as algebraic sums of these harmonics 
and frequencies u and 2u . 

As we approach the resonant field further, the shape of absorbed power peaks experiences 
rapid qualitative changes. Figure 56 shows the result of transformation of spike-like signals 
into typical saw-tooth temporal series for H = 8.47 kOe. In addition to the increase in the 
amplitude and relative height of frequency peaks, the emergence of linearly increasing and 
decreasing segments on the time dependence of absorbed power is also worth noting. It 
is remarkable that such oscillations are almost indistinguishable from classical relaxation 
oscillations that are frequently encountered in electrical engineering. 

A subsequent decrease in the field leads to the tendency to the formation of periodic 
rectangular pulses of absorbed power. Signals of such a shape are shown in Fig. 5c for 
H = 8.44 kOe. It should be noted that the amplitude of oscillations does not increase any 
longer, while the periodicity is enhanced, which is manifested in the frequency spectrum. 

Relaxation oscillations become completely chaotic for field values close to resonance. Fig- 
ure 5d shows the corresponding temporal series of absorbed power and a typical "grass-like" 
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frequency spectrum for H = 8.39 kOe. It can be seen that the amplitudes of oscillations are 
much smaller than those in Fig. 5c, and the frequency distribution of oscillations has become 
almost continuous with a sharp decrease in the maximum peak heights to the amplitude of 
the 50-Hz peak of the fundamental harmonic. 

As the field decreases further from the resonant value, relaxation oscillations again ac- 
quire the spike-like shape, being essentially nonlinear. Figure 5e shows for H = 8.34 kOe 
the temporal series for such anharmonic oscillations of absorbed power and their frequency 
spectrum with clearly manifested peaks of multiple harmonics. A distinguishing feature of 
these oscillations is that their fundamental frequency is almost half the frequency of similar 
spike-like oscillations presented in Fig. 5a. The frequency of relaxation oscillations in general 
decreases as the field decreases to the resonant value, and starts increasing after the passage 
of the resonance peak. 

In the magnetic field scanning in the opposite direction (i.e., upon its increase), the 
regimes described above appear in the reverse order, but a hysteresis loop takes place in 
complete accord with the picture shown in Fig. 1. 

It was mentioned above that the selection of other values of power (for example, the 
maximum power P = dB) followed by scanning in the static magnetic field results in 
chaotic nonlinear oscillations whose frequency structure contains higher harmonics of the 
fundamental frequency as well as subharmonics against the background of a high-intensity 
continuous noise spectrum (see Figs. 2a and 6). 

A comparison of temporal series also leads to the conclusion concerning clearly manifested 
temperature dependence of the degree of stochastization of oscillations. The higher the 
temperature, the higher the noise level in the oscillatory spectra and the extent of their 
nonregularity, and vice versa. At low temperatures, we could observe relaxation oscillations 
in the form of nearly rectangular pulses (such a mode was realized for P = 5 mW, H = 
8.3 kOe, and T = 1.7K). 

Another interesting feature is the observation of the regime of an abrupt and virtually 
complete disappearance of free oscillations with simultaneous doubling of the period of 
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nonlinear reference oscillations and the emergence of their subharmonic at the frequency 
25 Hz. We can try to explain the latter effect from the point of view of the theory of chaos 
control and the emergence of higher (multiple) resonances. However, we shall not consider 
this problem here and analyze the structure of chaotic attractors of relaxation oscillations. 

ANALYSIS OF EXPERIMENTAL RESULTS 

The method of a nonlinear analysis of experimental temporal series has been worked 
out intensely during the last decade and is described in detail in a number of reviews and 
monographs.B'&il We shall use this method which involves the determination of the linear 
autocorrelation function for temporal series, the determination of "time delay", the con- 
struction of phase portraits of attractors in the corresponding "time delay" coordinates, the 
construction of interspike intervals and their analysis, the computation of the correlation 
dimension of attractors, the determination of the noise contribution to temporal series, the 
source of the noise and possibilities of its reduction and the discussion of theoretical models 
of the observed chaotic oscillations. 

The temporal series is a discrete set of values of the physical quantity (the absorbed 
power V(t n ) in our case), measured in equal intervals of time. A traditional characteristic 
of temporal series of signals is the linear autocorrelation functional 

N 

jj £ [s(m + T) - s\[s(m) - s] 

Cl(t) = ^ , (2) 

£ E Hm) - s] 2 

m=l 

where the average value of the signal s(m) is defined in the standard manner: 

N 

s = n £ s(m). 

m=l 

Since we usually subtract the average value of series from the initial series in an analysis 
of spectra, we calculated autocorrelation function for time dependences presented in Fig. 3 
for modified series s — 0. As a function of r, it exhibits qualitatively identical behavior 
for all values of power: this is an oscillating function with a slowly decreasing amplitude. 
The period of these oscillations coincides with the fundamental period of oscillations of the 
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signal being measured. In an analysis of nonlinear signals, autocorrelation function is also 
useful for estimating "time delay". It is chosenlll'il equal to the value of r for which the 
autocorrelation function vanishes for the first time. In our measurements, this time delay is 
approximately equal to a quarter of the fundamental period of observed oscillations. 

Figure |6] shows the dependence of time r on the driving power P. (The unit of measure- 
ments of t is the principal interval At = 4.9 ms of our temporal series.) It can be seen that 
this dependence obviously correlates with the dependence of the frequency of oscillations on 
the driving power presented in Fig. 4 and confirms the existence of a threshold transition 
from one regime of chaotic oscillations to another. It should be noted that the period of the 
correlation function corresponds to oscillations of frequency 50 Hz for low powers and to the 
fundamental period of saw-tooth oscillations for the maximum power. 

The obtained value of r can now be used for plotting phase portraits of nonlinear os- 
cillations. For this purpose, we shift the temporal series by r and plot the dependence of 
V(t n + r) on V(t n ). These functions are just the time-delay coordinates. For the tempo- 
ral series corresponding to the maximum power in Fig. 3, the phase portrait is shown in 
Fig. [7|a (the value of r is chosen equal to 49 ms, and the average value of absorbed power is 
subtracted from the given series). It can be seen that the process is periodic on the whole 
and occurs in several stages with their own characteristic times. In order to obtain a more 
detailed concept of the attractor structure, we constructed a ID map from the sequence 
of minimum values of the Poincare sections of the given attractor. These values were de- 
termined as negative values of V(t m + r) taken at instants t m for which V(t m ) = 0. This 
dependence is shown in Fig. 7b and demonstrates the existence of the internal structure of 
the attractor and an obviously large contribution of noise. 

A detailed analysis of the time dependence V(t) indicates that the noise contribution is 
not additive. Indeed, irregular amplitude jumps as well as periodic oscillations of frequency 
50 Hz have different values at different stages of variation of the function V(t). This indicates 
a nonlinear enhancement of both factors and their participation in the chaotic process. 

In order to describe the chaotic behavior of relaxation oscillations and the effect of noise 
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on them quantitatively, we consider a sequence of time intervals between adjacent peaks of 
the signal and the sequence of maximum values of peak amplitudes as characteristics of this 
process. 

These dependences of amplitude peaks and interspike intervals for the series under inves- 
tigation are compared with the relevant sequences for the series shown in Fig. 5a. It can be 
seen that the amplitude peaks of oscillations of the absorbed power (Figs. [8|a and 6) behave 
quite chaotically in the vicinity of a resonance and far away from it. At the same time, the 
interspike intervals (Figs. 8c and d) exhibit a clearly manifested tendency to a quasiperiodic 
mode far away from the resonance (Fig. 8d) , but random forces acting on the system result 
in the stochastization of oscillations, which is accompanied by chaotic jumps in the period 
of oscillations between its four principal values. 

We can try to determine whether a noise is stochastic or dynamic by calculating the cor- 
relation dimension D of the attractor under investigation. This quantity is defined through 
the pair correlation integral: 



where N is the number of measurements, r the correlation radius, y m {i) 
the vector of dimensions m in the embedding space, whose coordinates are 
{V(ti)V(ti + t), ...V(ti + (m — 1)t)}, and 6{r) is the theta function. This function in fact 
determines the number of pairs of vectors in the m-dimensional space, the separation be- 
tween which is smaller than the preset distance r. While determining the distance, we 
presume that the cells into which the phase space is divided have the cubic shape.il The 
dimension D is the limit of the expression 



and is usually calculated on the interval r in which the values of the correlation function are 
not very small. 

The sequence of calculated correlation sums for the initial temporal series corresponding 
to the maximum power in Fig. 3 is presented in Fig. || on logarithmic scale (the curves 






(4) 
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correspond to the variation of m from 1 to 9 from top to bottom). Numerical differentia- 
tion reveals a flat segment according to which the dimension of a strange attractor can be 
estimated. It was found that it is slightly larger than two (2.25 ±0.1), but the strong effect 
of dynamic noise following from the characteristic increase in the steepness of the curves for 
ln(r) < 2.5 does not allow us to establish the existence of exact limit. We are inclined to 
interpret the latter quantity as the dimension of a regular attractor. On the other hand, 
the slope of the curves in the region 1 < ln(r) < 2.5 also demonstrates the tendency to 
a limit that can be estimated as 4.9 ± 0.1. Such a limiting value can be regarded as the 
total dimension of the attractor, containing the contribution from a regular attractor and 
a deterministic noise. A slight increase in the dimension for large values of m for small r 
is associated with the contribution of "white" instrumental noise. The above comparative 
analysis of functional dependences of temporal series and their spectra also confirms this con- 
clusion. Thus, the analysis of the correlation dimension leads to the conclusion concerning 
the dynamic nature of noise in the relaxation oscillations under investigation and indicates 
a multidimensional chaotic dynamics and, generally speaking, multimode excitations in the 
resonance system in question. The extent of its stochastization is quite high, which follows 
from the estimate of correlation dimension in the range of small r. It should be noted that 
recent investigations of parametric resonance in a related metallorganic antiferromagnet0 
also confirms the deterministic origin of noise in these compounds. 

The existing theories of relaxation oscillationslli'iilii make it possible to describe the 
emergence of high- dimensional chaos on the basis of multimode models. The main idea 
behind the mechanism of emergence of relaxation oscillations can be demonstrated even by 
using a two-mode model in which it is assumed that the resonance excitation conditions 
are satisfied for one mode and are not observed for the other mode. Under the action of 
pumping, such a pair of coupled nonlinear oscillators reproduces quantitatively the behav- 
ior of relaxation oscillations. A quantitative theory of this effect for antiferromagnets has 
not been developed as yet. However, a theoretical description of this phenomenon in the 
approximation of two spins simulating the sublattices subjected to resonant transverse and 
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longitudinal pumping appears as promising. Consequently, chaotic relaxation oscillations 
can be described qualitatively as the dynamics of a nonlinear oscillator under resonance 
conditions, but under the action of certain random forces (the inclusion of the effect of the 
second oscillator). Such a system may have at least two stable states the transition between 
which can lead to the emergence of spike-like and saw-tooth time dependences of absorbed 
power. The features and diversity of existing chaotic modes are obviously determined by the 
time of residence of the system in the equilibrium states and the rate of transient processes. 
Such a system can obviously have a high degree of stochastization, an attempt to create 
regular attractors in it will lead to regimes in which such attractors coexist with a well- 
developed dynamic noise. However, a quantitative theory of such chaotic oscillations should 
apparently be constructed on the basis of a multimode model according to the numerical 
analysis carried out by Moser et al.S 

The theorem on dimension that has been formulated recently for systems with a dynamic 
noise indicates in its simplified formulation the additivity of the dimension of a regular 
attractor and a noise.0 In this sense, these can be separated, and the question of elimination 
of noise from a signal, at least the noise reduction, and isolation of a regular signal from the 
data on temporal series appears as justified. 

There are effective methods of suppression of noise effect on the useful signal.il These 
methods are extremely effective for reduction of external additive noise, but their application 
in the case of a dynamic noise should be verified in each specific case. The algorithm of 
purification of a signal in the simplest form can be described as follows. In the chosen 
embedding space whose dimension is larger than the sum of the predicted dimension of the 
regular attractor and dynamic noise, the nearest neighbors of the preferred vector of state 
are selected, and its central coordinate is averaged over the values of relevant coordinates 
of the found neighbors. The obtained sequence of new data is the result of one iteration 
that can be repeated. Such an algorithm can be optimized as well as the choice of required 
parameters (correlation radius, dimension, etc.; see Ref. Such an algorithm will be used 
below for analyzing the chaotic temporal series measured by us. 
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On the other hand, the above algorithm in the simplest form includes the conventional 
method of data averaging over nearest and next to nearest neighbors in the series. In this 
case, the dimension of the embedding space is equal to unity, and the number of neighbors is 
fixed. It can easily be verified that in spite of its very simple form, the procedure operates as 
a high-frequency filter and does not change the complex low-frequency spectrum of chaotic 
oscillations. We shall apply this procedure also to analyze the results. 

We chose the object of investigation in the form of a chaotic attractor obtained from the 
temporal series presented in Fig. 26. The results of analysis are shown in Fig. |10| (it should 
be noted that the average value is subtracted from the terms of the series). Figures 10a and b 
show the phase portrait on the plane {V(t), V(t+r)}, where r = 29.5 ms, and simultaneously 
the map for the data on the Poincare cross section (see above), while Figs. 10c, d, e, and 
/ contain the dependences for the data "corrected" by the method of averaging and the 
optimized method of noise reduction described above respectively (averaging was carried 
out twice over five points, and tenfold iterations were used in the optimized method). 

Figures lOe and / clarify the internal structure of a regular attractor. After the effect 
of noise becomes weaker, its phase portrait resembles a strange multiband attractor. The 
analysis of correlation dimension makes it possible to characterize quantitatively both the 
regular attractor as well as the residual contribution of deterministic noise. The results of 



analysis are presented in Fig. [IT]. It should be noted that prior to calculation correlation 
functions for the temporal series under consideration, we initially normalized all values to 
a unit interval by the formula V(ti) = (V(ti) — V m i n )/{V max — V min ), where V min and V max 
are the minimum and maximum values of the signal in the series. It was found that the 
dimensionality of a regular attractor can be estimated as 2.15 ± 0.05, and the total dimen- 
sionality with the contribution of deterministic noise as 3.25 ± 0.05. Spectral analysis of 
these "improved" results also indicate that the quantitative contribution of noise remained 
quite large, and the resultant attractor possesses a high dimensionality as before. 

Thus, chaotic dynamics in the nonlinear antiferromagnetic resonance in low-dimensional 
antiferromagnets is high-dimensional, the extent of stochastization of oscillations is high, 



and noise has a deterministic origin and serves as a decisive factor in nonlinear dynamics of 
these magnets. 

Finally, we formulate the following conclusions following from our analysis. 

1. Peculiarities of a transition to chaos by "irregular periods" in a 2D metallorganic 
antiferromagnet with an "easy axis" type anisotropy are experimentally observed and 
studied in detail under conditions of nonlinear antiferromagnetic resonance. 

2. It is shown that relaxation oscillations of absorbed power are generated for very low 
energy levels of microwave field and have a low frequency of fundamental harmonic 
(of the order of a few hertz). No multiple harmonics are observed experimentally at 
kilohertz and higher frequencies. 

3. Relaxation oscillations at low values of driving power exist in the form of generally 
periodic sequence of spike-like peaks of absorbed power. The frequency spectrum 
contains components of fundamental frequency corresponding to the emergence of 
spikes as well as multiple harmonics, which demonstrates a nonlinear nature of the 
process. 

4. As the pumping amplitude increases, the phenomenon of period doubling is observed 
in the time dependence of absorbed energy of microwave field. The shape of the 
signal is simultaneously transformed from the spike-like to the saw-tooth type having 
segments with linearly increasing and linearly decreasing absorption. An analysis of 
the frequency-amplitude dependence of oscillations and their linear autocorrelation 
function gives quantitative characteristics of this transition. 

5. A similar effect is observed at a fixed level of pumping, but upon a change in the value 
of static magnetic field near its resonant value. 

6. With increasing power, relaxation oscillations become chaotic. The spectrum of such 
oscillations is continuous and has a "grass-like" form, but the peaks of fundamental 
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harmonics are still distinguishable. The phase portrait of these oscillations has the 
form of a strange attractor experiencing a strong influence of noise. Stochastization of 
oscillations, however, is not a result of influence of an additive instrumental noise. 

7. The quantitative characteristics of such a strange attractor are calculated. The one- 
dimensional map corresponding to the given attractor demonstrates a tendency to 
regular movement in spite of chaotic time dependence of relaxation oscillations. An 
analysis of correlation dimension indicates the high-dimensional chaos dynamics and 
the deterministic nature of noise in the magnetic system under investigation. The 
possibility of formal separation of the regular movement and the noise contribution 
with the help of nonlinear methods of noise reduction being developed is considered. 

8. The applicability of the theoretical model of finite number of coupled spins under the 
action of the parametric and transverse pumping to the construction of a quantitative 
theory of the scenario of transition to chaos by "irregular periods" is discussed briefly. 
The transition can be regarded as a universal phenomenon in low-dimensional ferro- 
and antiferromagnets under nonlinear resonance conditions. 
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FIGURES 

FIG. 1. Amplitude-field dependences of nonlinear antiferromagnetic resonance for a 2C4Mn 
crystal for different values of driving power. The bold segment on the upper curve in the vicinity 
of the high-field peak denotes the range of relaxation oscillations. The hysteresis loop observed for 
a power exceeding — 5dB increases with driving power. 

FIG. 2. Chaotic oscillations of absorbed power and their frequency spectrum for a driving 
power of 5 mW for a resonant curve with closely spaced peaks in fields H m = 8.37kOe(a) and 
8.4 kOe (6). 

FIG. 3. Evolution of time dependences of absorbed power for a change in the microwave power 
level from —10 to dB for a static magnetic field H = 8.4 kOe. 

FIG. 4. Frequency of relaxation oscillations as a function of driving power. Solid and open 
circles correspond to an increase and decrease in the driving power respectively. The threshold 
effect of oscillation period doubling is of the hysteresis type. 

FIG. 5. Temporal series of absorbed power and their spectra for different values of static mag- 
netic field for the driving power P = — ldB : (a) H = 8.51 kOe, relaxation oscillations with a 
spike-like structure and a high stochastization level; (b) H = 8.47 kOe, the result of transformation 
of spike-like signals into saw-tooth signals; (c) H = 8.44 kOe, saw-tooth relaxation oscillations with 
linearly increasing and decreasing segments; (d) H = 8.39 kOe, chaotic temporal series with an 
intense "grass-like" frequency spectrum; (e) H = 8.34kOe(e), nearly regular anharmonic oscilla- 
tions of absorbed power of a spike-like shape. All the spectra contain the peak of the fundamental 
frequency of relaxation oscillations of the order of several hertz and its higher harmonics as well 
as the peak of the reference frequency 50 Hz and combination frequencies. 

FIG. 6. Dependence of "time delay" r on the pumping power P. The value of r is measured 
in units of the temporal series period At = 4.9 ms. 
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FIG. 7. Dynamic characteristics of relaxation oscillations for the driving power P = OdB and 
H = 8.4 kOe, whose time dependence is shown in Fig. 3: (a) phase portrait in the time-delay 
coordinates (r = 49 ms) and (b) one- dimensional map constructed from the sequence of minimum 
values of the Poincare cross section of the phase portrait taken at the instants of time t m at which 
V(t rn ) = 0. 

FIG. 8. Amplitude peaks and interspike intervals for relaxation oscillations for driving power 
P = OdB and H = 8.4 kOe (see Fig. 3) (a,c), and P = ldB and H = 8.51 kOe (see Fig. 5a) (b,d). 

FIG. 9. Dependences of logarithms of correlation integrals on the logarithm of the distance 
between vectors in the m-dimensional embedding space calculated from the initial integral data 
for the attractor presented in Fig. 7. 

FIG. 10. Phase portraits and one-dimensional map for relaxation oscillations for the driving 
power 5 mW in the field H = 8.4 kOe in the case of closely spaced peaks on the resonant curves (see 
Fig. 26) for original data on temporal series (after the subtraction of the average value) (a, 6), for 
doubly averaged data for five nearest neighbors (c, d), and for data "improved" by the optimized 
method of noise suppression as a result of ten iterations (e, /). 

FIG. 11. Dependences of logarithms of correlation integrals on the logarithm of the distance 
between vectors in the m-dimensional embedding space calculated for data "improved" by the 
optimized method (see Figs. lOe, /). 
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